Nonlocal Kinetic Equation and Simulations of Heavy Ion Reactions 
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A kinetic equation which combines the quasiparticle drift of Landau's equation 
with a dissipation governed by a nonlocal and noninstantancous scattering integral 
in the spirit of Enskog corrections is discussed. Numerical values of the off-shell 
contribution to the Wigner distribution, of the collision duration and of the collision 
nonlocality are presented for different realistic potentials. On preliminary results 
we show that simulations of quantum molecular dynamics extended by the nonlocal 
treatment of collisions leads to a broader proton distribution bringing the theoretical 
spectra closer towards the experimental values than the local approach. 

Introduction 



OCPne of the long standing problems in nuclear physics is to find the equation of state of nuclear matter 

In the absence of any direct measurement, it is hoped that the equation of state can be deduced 

'from heavy ion reactions via dynamical simulation of the fragmentation scenario. Most simulations 
<-> : 

^•ely on the local and instantaneous treatment of binary collisions as they appear in the Boltzmann 
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The arguments of distributions / and energies e are shortened as f\ = f a (k,r,t), fi = fb(p,r,t), 
f 3 = f a {k — q,r,t), and / 4 = fb(p + 1, r ,t), with momenta k,p,q, coordinate r, time t, and spin and 
isospin a, b. The local picture of the collision is reflected by the same coordinate r at all distributions, 
the instantaneous by the same time t. 

A real binary collision is neither local nor instantaneous. A nonlocal picture of a collision is 
schematically drawn in figure |l[ Let us introduce individual nonlocal corrections step by step for 
simple but useful models. The first nonlocal corrections have been introduced by Enskog for the 
classical gas of hard spheres 0. The collision of hard spheres is instantaneous so that each trajectory 
is broken only at a single point. Accordingly, in figure A 3 = and A 4 — A 2 = 0, i.e., Af = and 
Aj, = 0. At the instant of the collision, the particles are displaced by the sum of their radii in the 
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Figure 1: A nonlocal binary 
collision. 

direction of the transferred momentum what is described by a nonzero vector Ahs- Corresponding 
changes in the kinetic equation enter the position of the ongoing particle, f 2 = fb(p, r — Ahs,^) and 
f& = fb(p + q,r — Ahs,£) while other arguments remain unchanged. The equation of state evaluated 
from the kinetic equation with the nonlocal scattering integral is of the van der Waals type covering the 
excluded volume 0, [J. For nuclear matter, Enskog's corrections have been first discussed by Malfliet 
(U and recently implemented by Kortemayer, Daffin and Bauer ||. 

Another simple model one obtains assuming that colliding particles form an unstable molecule 
of an average time of life A t called the collision duration or the collision delay. Neglecting the size of 
the molecule, A 2 = and A 4 — A 3 = 0, i.e., Ahs = and = 0, the only nonzero displacement 
Af measures a distance between points where molecule forms and breaks up. This distance is given 
by Af = A t f mo1 with the molecular velocity t> mo1 = (k + p)/(m a + nit,). Corresponding changes in the 
kinetic equation enter the time argument and positions of final states, f% = f a (k — q,r — A{,t — A t ) 
and / 4 = fb(p + q, r — A f , t — A t ). In the equation of state, the collision duration results in the same 
kind of terms as the presence of stable molecules. The finite duration of nucleon-nucleon collisions 
and its thermodynamic consequences has been studied by Schmidt, Ropke and Schulz ||, its effect 
on the pressure has been discussed only recently by Danielewicz and Pratt J7|. The noninstantaneous 
scattering integral and its consequences for the linear response has been first discussed for electrons in 
semiconductors scattered by resonant levels |§. 

In a real collision, the two particles keep a finite distance, Ahs 7^ 0, and interact for a finite 
time, At ^ and Af ^ 0. Moreover, particles rotate one against the other as described by the 
rotation displacement A^. Consequently we obtain the following scenario of Fig. P Two particles 
approach until they reach a distance A 2 . Then they form a molecule living for At and traveling over a 
distance Af . During this propagation the molecule rotate as given by A^. Collecting all three shifts we 
obtain the nonlocal and noninstantaneous kinetic equation (|6]) derived in || [H| with the help of the 
method introduced in ||. The resulting arguments of the kinetic equation read finally f\ = f a (k,r,t), 




/a = fb(p,r- A 2 ,t), / 3 = / (fc - g,r - A 3 ,t- A t ) and / 4 = fb{p + q,r - A 4 ,t - A t ). 

While the above microscopic picture of nonlocal and noninstantaneous isolated collisions is intu- 
itively clear, it is less transparent how to define the same corrections for quasiparticles which carry a 
part of the interaction in the quasiparticle reconstruction of their energies and wave function norms. 
This question requires a systematic approach as it was first presented in M. This derivation follows 



Baerwinkel [11] in starting from nonequilibrium Green's functions and keeping all gradient contribu- 
tions to the scattering integral, but instead of the quasiparticle approximation, the extended quasipar- 
ticle approximation is used. Here we use numerical results to discuss the key steps and consequences 
of this approach. 



2 Kinetic equation 

We start our derivation of the kinetic equation from the quantum transport equation for the nonequi- 



librium Green's functions first obtained by Kadanoff and Baym, see [12 
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G 1 - ReS, G<] - i [ReG, £<] = \ {G>, £<} - \ {G<, £>} , (2) 



where [, ] and {, } denote commutators and anticommutators, ReG = \{G R + G A ) is the hermi- 
tian part of the propagator. The center of interest is the particle correlation function G < (1,2) = 
TR (ptft(2)tf(l)). Its time evolution, however, requires to know the accessible states given by the 
hole correlation function G > (1, 2) = TR (pty (1)^(2)), and a dynamics of interactions specified by the 
selfenergy S. Individual terms in (0) have specific physical content. The G$ 1 describes a free motion 
of particles and is renormalized by ReS. The ReG describes the off-shell motion after the collision. 
The first and second anticommutators are the scattering- in and -out. 

The dynamics of interaction reflects selected models and approximations. For simplicity we 
assume that protons and neutrons are of equal mass m, interact via an instant potential V, and there 
is no spin-flipping mechanism. As common, the self-energy is constructed from retarded and advanced 
two-particle T-matrices T R ' A in the Bethe-Goldstone approximation \T2, |TB| as 

£<(1, 2) = T R {13, 56)T A (78, 24)G>(4, 3)G<(5, 7)G<(6, 8), (3) 

and E > is obtained from @ by an interchange >^<. Numbers are cumulative variables, 1 = 
(ti,ri,ai), and bars denote internal variables that are integrated over. Missing commas in argu- 
ments signal that the time arguments are identical, e.g., £g = t% and tg = ig, the T-matrices are 
thus double-time functions. The approximations of set are specified in the mixed representation, 
[off-diagonal elements in spin and isospin are excluded, a\ = a 2 = a] 

G<(1 ' 2) = / ^(^a^^" 2 ^^^" 2 ^" i^r,t) r= ^ t=h ^ . (4) 



2.1 Off-shell motion 



Equations (|2|-[3|) completed with the Dyson equation for G R and the ladder equation for T R form a 
closed set for G < . This set is converted to equations for the quasiparticle distribution / with the help 
of the extended quasiparticle approximation [|14], [L3|, [L5|, [L6[ 

(5) 
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where G\ )LU = G a (u>, k,r,t) and similarly X. The first term is singular and provides the dominant 
quasiparticle contribution on the energy shell. The second term is regular and contributes out of 
the energy shell. The approximative form of this off-shell contribution is consistent with the lowest 



approximation of the wave-function renormalization, z\ 
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Figure 2: The equilibrium occupation of mo- 
mentum states given by Wigner's distribu- 
tion p (full line) is compared with the quasi- 
particle distribution / (long dashed line) and 
with the extended quasiparticle approxima- 
tion (short dashed line) at the temperature 
T = 10 MeV. 



extended quasiparticle approximation (|S|) to convert the transport equation (|2|) into a kinetic equation, 
it is useful demonstrate its physical content. The amplitude of the off-shell contributions and the need 
to treat them separately as in (|5|) can be seen in figure |2| where we compare the Wigner distribution 
Pi = 2~ J dooGf^ obtained from the exact equilibrium correlation functions G < with its extended 
quasiparticle approximation (||) and the plain quasiparticle approximation Gf w = /i27r5(a; — 6\) which 
yields the Fermi-Dirac distribution. We note that the simple Yamaguchi potential []l7| has been used 
in the T-matrix for this comparison. One can see that the off-shell contribution given by the difference 
between the Wigner and the Fermi-Dirac distributions is not small, in particular at the high momenta 
region where the power-law off-shell tails always dominate over the exponentially falling quasiparticle 
distribution. Formula (|5|) provides inevitable and sufficiently precise off-shell corrections. 

Out of equilibrium one has to evaluate S < which is similar to an evaluation of the scattering 
integral. From the extended quasiparticle approximation one then obtains the high-momenta tails 
found by Danielewicz |18[ and Kohler ]19|] in direct numerical treatments of non-equilibrium Green's 
functions. 



2.2 Quasiclassical trajectories 



Now we substitute (|5]) into (0) and neglect all gradient terms but linear. Within nuclear physics, the 
gradient expansion of the self-energy (§) is commonly omitted since it is a tedious task. It results in 
one nongradient and nineteen gradient terms that are analogous to those found within the chemical 
physics | f20| , |2~T|j . All these terms can be recollected into a nonlocal and noninstantaneous scattering 
integral that has an intuitively appealing structure of a nonlocal Boltzmann equation (fj) 
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with Enskog-type shifts of arguments [j9j |10|: fx = f a (k, r, t), f 2 = f b (p, r-A 2 , t), f 3 = f a {k-q-A K , r- 
As,t — A t ), and f^ = fb(p+Q — A^, r — A4, t — A t ). In agreement with |20|, all gradient corrections 
result proportional to derivatives of the scattering phase shift cf) = Im lnT^(f2, k,p, q, t, r), 
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and A r = \{A 2 + A 3 + A 4 ). After derivatives, A's are evaluated at the energy shell Vt — > e 3 + £ 4 . 

For the purpose of discussion, it is advantageous to link the quantum displacements to 
intuitively more appealing hard-sphere and rotation shifts by relations obvious from figure [l] 



A HS = -(A 4 - A 3 + A 2 ), 



-(A 4 
2 V 



A, - A 2 



(8) 



For the collision of two isolated nucleons, it is possible to show that Ahs points in the direction of the 
transferred momentum q. Similarly it follows that the rotation shift is orthogonal to Ahs and stays in 
the collision plane. 

In figure ^ we plot the delay time and the amplitude of the hard-sphere shift for different deflection 
angles versus lab energy. The T-matrix is evaluated with different potentials, Bonn (A-C) f22|, Paris 
[p3| and separable Paris f2l |, concerning partial wave coupling up to D- waves [^5]. The forward angle 
delay time has a negative minimum at small energies indicating an attractive behavior. For very small 
energies the delay time is rapidly decreasing to high negative values reflecting the occurrence of weakly 
bound states. The sharp jump for 6 = 90° is caused by a resonant character of this scattering channel. 
Its value is misplaced and exaggerated within the separable approximation. The hard core shift Ahs 
behaves regularly at low energies. The resonant scattering at 9 = 90° appears as the increase of the 
amplitude. Again, the separable potential exaggerates its value. Characteristic values A t ~ 1 fm/c 
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Figure 3: The collision delay A t (left) and the hard core shift |Ahs| (right) for isolated nucleons in the 
barycentric coordinate system for different deflection angles 9 and interaction potentials. 

and Ahs ~ 1 fm show that the noninstantaneous and nonlocal treatment of binary collisions can be 
important in heavy ion reactions. 



3 Implementation in heavy ion codes 

The selfconsistent evaluation of all A's would be as demanding as the full Green's function treatment 
of the system. We employ two kinds of additional approximations. First, following the approximations 
used within the Boltzmann equation, we neglect the medium effect on binary collision, i.e., use the 
well known free-space T-matrix. Second, we rearrange the scattering integral into an instantaneous 
but nonlocal form. To eliminate the time delay from distributions we neglect collisions on the time 
scale A t and shift arguments of distributions along the trajectory of a particle. In the differential form 
this step derives as 

f a (k-q-A k ,r-A 3 r ,t-A t ) = f a (k-q-A k ,r-A 3 r ,t)-^-A t 

= fa(k-q-A k -^A t ,r-A 3 + ^A t ,t\. (9) 

In the last step we have used the collision-free kinetic equation ^ + — ^r^f — 0. With 

approximation ([]) and similar for / 4 , the scattering integral (|6|) corresponds to an instantaneous event 
at time t. The collision remains nonlocal with modified space displacements of particles at initial states 

A 3 ,4 = A 3 ,4 - -^A*. (10) 



The instantaneous approximation brings further simplifications following from conservation laws. 
During the instantaneous process, mean fields have no time to pass any momentum and energy to the 
colliding pair. Indeed, assuming the effect on colliding particles only via mean fields, from ([?[) one finds 
A& = — ^^-A t so that the momentum gain vanishes in @. Similarly, the energy gain Ag vanishes 
when arguments of quasiparticle energies in the energy conserving 5 function are brought to the same 
time instant. Finally, in agreement with the continuity of the center of mass motion, one finds that 
A2 = A3 + A4. The scattering- in thus simplifies as 

E / J^J^3 2n5 ^ + £ 2~e 3 - 5 4 ) \T ab \ 2 (k - p, q) hh{\ - h){\ - f 2 ), (11) 

where new arguments of energies and distributions are shifted only in space, f\ = f a (k,r,t), f 2 = 
fb(p, r — A 2 , r, t), f 3 = f a (k — q, r — A 3 ,t) and = fb(p + q, r — A 4 ,t). The scattering-out is similar. In 
the T-matrix we have reduced arguments to those which are relevant for free-space collisions. Finally, 
we would like to stress that the amplitudes of displacements are not fitting parameters or a matter of 
an educated guess but evaluated from the T-matrix, i.e., from an interaction potential. 
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Figure 4: The exclusive QMD proton spec- 
tra for central collision of 129 Xe^ 119 Sn at 
50 MeV/A with and without non-local cor- 
rections. The data are extracted from recent 
INDRA experiments. The non-local correc- 
tions bring the spectrum towards the exper- 
imental values. 
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In order to investigate the effect of nonlocal collisions in realistic situations, we have evaluated 
A's from the separable Paris potential and implemented the nonlocal scattering integral of ( |Tl~D in 
a QMD code for the central collision of 129 Xe— > 119 Sn at 50 MeV/A. Figure |] shows the exclusive proton 
spectra subtracting the protons bound in clusters. This procedure is performed within a spanning tree 
model which is known to describe a production of light charged cluster in a reasonable agreement with 
the experimental data. Within the local approximation, however, the remaining distribution of high- 
energy protons is too low to meet the experimental values. As one can see, the inclusion of nonlocal 



collisions corrects this shortage of the QMD simulation. The increase in the high-energy part follows 
from an enhancement of a number of collisions at the pre-equilibrium stage of the reaction while later 
stages are not strongly affected ||26| . Accordingly, the production of light clusters is rather insensitive 
to the nonlocal corrections. The improvement of the proton production is thus not on cost of worse 
results in other spectra. 



4 Thermodynamic properties 

The production of high-energy particles documents a vital role of nonlocal treatment far from equilib- 
rium. Their role can be best seen on thermodynamic observables like density n a of particles a, density 
of energy S, and the stress tensor Jij which conserve within the nonlocal and noninstantaneous kinetic 
equation (|g). Integrating @ over momentum k with factors e±, k and unity one finds || [T(| that each 
observable has the standard quasiparticle part following from the drift 

£ "' = ?/(lp^ + ^/||^ 1 + ,^, P ,o)/ 1 / 2 , 



a, ft 

qp 



and the A-contribution following from the scattering integral 

= I E/ Kp+?) A* + (k-q)iA 3j - Pi A 2j ] , 

where P = \T a b\ 2 2ir5(ei+E2 — £3—£4)fif2(l- — f3 — f4)- The arguments denoted by numerical subscripts 
are identical to those used in ([l]), for all A's are explicit. 

The density of energy £ = S qp + AS alternatively results from Kadanoff and Baym formula, 
£ — J2 a I J0js I fjff (w + ^n) ^ai^i r i Oj w ^h G < in the extended quasiparticle approximation (||). 
Its complicated form, however, shows that 8 cannot be easily inferred from an eventual experimental 
fit of the kinetic equation as it has been attempted in . The conservation of £ generalizes the result 
of Bornath, Kremp, Kraeft and Schlanges [HJ restricted to non-degenerated systems. The particle 
density n a = n qp + An a is also obtained from (|]) via the definition, n a — J 2 ^10f ( ^ < - This confirms 
that the extended quasiparticle approximation is thermodynamically consistent with the nonlocal and 
noninstantaneous corrections to the scattering integral. 

For equilibrium distributions, formulas (|T2|) and ([13]) provide equations of state. Two known 
cases are worth to compare. First, the particle density n a = n qp + An a is identical to the quantum 



Beth-Uhlenbeck equation of state || [13|, [16[], where n^ p is called the free density and An a the correlated 
density. Second, the virial correction to the stress tensor has a form of the collision flux contribution 
known in the theory of moderately dense gases [Q, [|. 



5 Summary 

In this paper we have discussed the kinetic equation which is consistent with thermodynamic observ- 
ables up to the second order virial coefficient. This theory extends the theory of quantum gases and 
non-ideal plasma to degenerated systems. The amplitude of the contribution of the off-shell motion is 
demonstrated on the Wigner distribution which also shows the precision of the extended quasiparticle 
approximation. The contribution of the nonlocal corrections to the scattering integral is documented 
on a realistic study of the heavy ion collision within the quantum molecular dynamics. It should be 
noted, however, that the separable Paris potential used in the presented preliminary results is not fully 
reliable as seen from comparison of the collision delay and the hard core displacement. Nevertheless, 
this preliminary study shows that the nonlocal corrections can be evaluated from the T-matrix and 
incorporated into existing Monte Carlo simulation codes. 
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